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Abstract 

We examine the lepton-flavor violation caused by a Yukawa coupling matrix 
yu,ij for right-handed neutrinos in the supersymmetric standard model. We stress 
that decay rates for t ^ and /i — > 67 may reach the range to be accessible to 
near future experiments if left-right mixing terms in the slepton mass matrix are 
substantially large. 



A small, but non-vanishing neutrino mass, if any, is regarded as an important indica- 
tion of a new physics beyond the standard model. The most interesting candidate is the 
well-known seesaw model that explains very naturally the small mass for neutrino in 
terms of a large Majorana mass for a right-handed neutrino A^. This model has attracted 
many authors not only because of its natural structure but also because the presence of 
N would illuminate some of the deep questions in particle physics. 

From the phenomenological point of view, on the other hand, introduction of three 
families of the right-handed neutrinos A^* (where i and j are flavor indices) brings two 
new ingredients to the standard model; one is a new scale of the Majorana masses for 
A^*, Mj, and the other a new matrix for Yukawa coupling constants of A^*, Vu^ij- Thus, 
we have two independent Yukawa matrices in the lepton sector as in the quark sector. 
In general, a simultaneous diagonalization of the both matrices is quite accidental and 
then the addition of the new Yukawa coupling lepton-fiavor violation. In the 

standard model, however, the amplitudes for the lepton-fiavor violating processes at low 
energies are suppressed by an inverse power of Mj at least and hence we do not expect 
sizable rates for such processes as far as Mi is very large. The neutrino oscillation is a 
famous exception in this point of view. 

It has been, however, noted by several authors that in the supersymmetric standard 
model (SSM) the Yukawa coupling y^^ij of A^* generates off-diagonal entries in the mass 
matrices for sleptons through the renormalization effects, which leads to unsuppressed 
lepton-fiavor violations such as r ^ yU7, yU 67, and so on.Q The predicted rates for 
the individual processes depend on the unknown Yukawa matrix Uy/ij- If off-diagonal 
elements of yy^ij are very small like in the quark sector, however, the reaction rates for 
these processes are predicted too small to be accessible to the near future experiments. 

The purpose of this letter is to point out that left-right mixing terms in the slepton 

mass matrix may give larger contributions to the lepton-fiavor violating processes, which 

have not been considered in the previous literatures . We show, as a consequence, that 

the decay rates for t fi'j and fi ^ e'-f can reach indeed the range close to the present 

experimental upper limits Br(r — >■ n'-f) < 4.2 x 10~^ and Br(/i — 67) < 4.9 x 10~^^ ^ even 

^It has been pointed out very clearly that the similar lepton-fiavor violation occurs in the SUSY 
grand unified theories (GUTs). 
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when the off-diagonal elements of yu,ij are small as yv,2zlyv,z?, ~ 0.04 and yv^nlVv^zz ~ 0.01. 

In the SSM with three families of right-handed neutrinos A^* (i = 1 - 3), the superpo- 
tential W for the lepton sector is given by, 



(1) 



where L*, and A^* are the chiral multiplets for left-handed lepton doublets, right-handed 
lepton singlets and right-handed neutrinos, and Hi and H2 those for the Higgs doublets. 
For simplicity, we assume the mass matrix M/ jj is proportional to the unit matrix as 



M,,i,- = Mi 



( 1 



V 



(2) 



The analysis with the general case M/ = Mi^ibij is straightforward, but as long as 
Mii ~ Mij the conclusion will not be much different from that in the simplest case (^. 
We choose a basis on which y^^ij is diagonal at the electroweak scale m^. 



?/e,jj(/i = mz) = ye,Aj- 



(3) 



To demonstrate our main point, we assume tentatively that the Yukawa couplings of 
are identical to those of the up-type quarks, similar to the case of SO(IO) unification,^ 



Vu, 



(4) 



We suppose that these relations hold at the gravitational scale /i = Mq ~ 2.4 x lO^^GeV. 
With the basis where the Yukawa couplings for the down-type quarks are flavor-diagonal 
at the electroweak scale, yd,ij = yd,iSij, the yu^ij is written as 



KM 



/ Vu 



V 



Vc 



\ 



yt J 



KM, 



(5) 



where is the Kobayashi-Maskawa matrix. The meaning of the Yukawa couplings of 
the up-type quarks yc and yt is the obvious one. 



^We do not use the other SO(10)-hke relation ye,ij = yd,ij for the down-type quarks, since it gives 
wrong results on the masses for the first and second families. 
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The seesaw mechanism induces small neutrino masses as0 



Mi 



( yu 






^ Vu \ 


Vc 






Vc Vkm 


V 


yt ) 




\ yt J 



(6) 



We see that the mixing matrix appearing in the neutrino oscillation is approximately 
identical with Vkm, due to the large hierarchy Vu ^ Vc ^ Vt 0- For a given M/, 
we determine the mass eigenvalues and the mixing angles for the neutrinos. We adjust 
M/ ~ lO^^GeV so that the mass for is predicted to be rrij,^ = lOeV,^ which lies in an 
interesting range in cosmology 0. To estimate the Yukawa matrix yu^ij{fi = Mg) we use 



the observed values m„(lGeV) = 4.5MeV, mdlGeV) = 1.27GeV l, mi = 174GeV ^ 
and Vkm- In our numerical calculations, we use the central value of each matrix element 
VKM,ii given in ref.[|]. 

We are now in a position to discuss the lepton-flavor violation. In the minimal SUSY 
standard model (MSSM), soft SUSY-breaking mass terms for sleptons have the general 
form, 

Csoit = -ml ,,f ^ - ml - (a,,,^'uH^ + h.c\ . (7) 



The lepton-fiavor conservation is easily violated by taking non-vanishing off-diagonal el- 
ements of each matrices and the sizes of such elements are strongly constrained from 



experiments. In the SUSY standard model based on the supergravity [0, it is therefore 
assumed that the mass matrices ml and ml are proportional to the unit matrix, and 
Ae^ij is proportional to the Yukawa matrix ye,ij- With these soft terms, the lepton-flavor 
number is conserved exactly. 



■^This equation holds only at tree-level. In our numerial analysis, we take into account the renormal- 
ization effects. 

^In this case, the mass for is predicted as ~ 10~'^eV, which is the right value for the MSW solution 
to the solar neutrino problem 0. However, the predicted mixing angle Ov^^^ ~ ^Cabibbo is much larger 
than that in the small angle solution of MSW ||]. Therefore, to account for the solar neutrino deficit, 
we must use smaller value for y^^i2- However, our main conclusions for r jj.'y and /j, e'y given in the 
text are unchanged, since the Yukawa coupling ?/j^i2 is almost irrelevant for these processes as noted in 
ref. i. 
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It is, however, not true if the effects of the right-handed neutrinos are taken into 
account [0]. The Yukawa couphng of neutrino y^^ijN^L^ H2 in eq.(|ID and the soft SUSY- 
breaking terms such as 

Aoft,. = -m%^^WN^^ -[A,,,,N'L^H2 + h.c.) (8) 

induce off-diagonal elements of through the radiative corrections. The renormalization 
group equation (RGE) for mj^ is given by 

drn?r [ drn?T Ir.oi, ,i o, ,4., o 



Hyl^lyu)^3 + '^{AlA,)ij], (9) 



where ((im| jj/(iln/i)MSSM represents the /3-function in the MSSM,^ and m|^2 is the soft 
SUSY breaking mass for the Higgs doublet if 2- As one can see easily in eq.(||), off-diagonal 
elements of m\ are induced by the renormalization effects if non- vanishing off-diagonal 
elements of and Ay exist. 

In order to obtain the slepton mass matrix m\ at the electroweak scale, we solve the 
RGEs for the full relevant parameters numerically. At the energy scale Mj < yU < Mq 
we use the RGEs derived from the MSSM with the right-handed neutrinos, and below 
the energy scale Mj the MSSM-RGEs without the right-handed neutrinos. For the soft 
SUSY-breaking parameters, we assume the following boundary conditions suggested in 



the minimal supergravity |Tl], 113 



mj^ijdj. = Mg) = mlSij, (10) 
Af^ij{n = Mg) = a mo yj^ij, (11) 

where iTijij is the soft SUSY-breaking mass matrix for sfermion / (with f = u, d, e and 
u), Af^ij the so-called A-parameter, y^ jj the Yukawa coupling constants for the fermion 
/, mo the universal SUSY-breaking mass, and a is a free parameter of 0(1). For gaugino 
masses, we use tentatively the relation implied by the SUSY GUTs ||13[| 



mG3 _ mG2 _ 3 tugi , . 

2 ~ 2 ~ r ^ 2 1 V / 

gi 92 5 gi 



^See ref. [n2[ for explicit formulae. 
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where mcs, '"^02, fnci, and (?3, g2, gi are the gaugino masses and gauge couphng constants 
for the gauge groups SU(3)c, 811(2)^ and U(l)y, respectively. 

The soft SUSY breaking mass matrix m| is determined by solving the coupled 
RGEs for all relevant parameters. Since tan/3(= {H2) / {Hi)) dependence of the result is 
relatively mild, we take tan /5 = 3 for the time being. The result for the larger tan (3 will 
be given later for comparison. For the case < 17102 < rno, the mass matrix w^i at the 
electroweak scale is given by 



mi(/i = mz) 



I (1.0-2.0) -(0.9 - 1.1) X 10"^ -(2.3 - 2.7) X 10-3 

-(0.9 - 1.1) X 10-^ (1.0-2.0) -(1.0 - 1.2) X 10-2 

V -(2.3-2.7) X 10-3 -(1.0- 1.2) X 10-2 (0.75 - 1.69) 

xm2, (13) 



where the magnitude of each entry becomes larger as mG2 larger. However, for the case 
^G2(Ai = ^z) > fno, ratios of the off-diagonal elements to the diagonal ones become 
smaller than those in eq.(|l3D, which give a more suppression of the lepton-flavor violation 
clS db result. 

We are now ready to calculate the reaction rates for various lepton-flavor violating 
processes. Let us first consider the diagrams considered in ref. [0] contributing to the 
process —>■ Ij'j. In this paper, we ignore the mixings in the neutralino and chargino 
sectors and consider that the bino and the winos are mass eigenstates. This approximation 
is justified if the SUSY-invariant higgsino mass is large, which is nothing but a situation 
we will consider in the present paper. 

From the diagrams in fig. 0, we obtain the amplitude for this process as 

Mih ^ Ijl) = ClJjPrU. AU. (14) 

where q and e are the momentum and polarization vector of the emitted photon. Here, 
the coefficient Cll is given by 



4 ■■"eL> 



+ —T^9lSe{ml2/mlj + -^glS^{mlJmlj L (15) 
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where m;., mg^ and m^;^ represent masses for the lepton charged cl and neutral sleptons 
ul, respectively,^ and the functions Se and are given by 

Se(x) = ——^ — -r{17x^-9x^-9x + l-6x^(x + 3)lnx}, (16) 
48(x — 1)^ 

Sy{x) = - ^ {-x^-9x^ + 93: + l + 6x(x + l)lna:}. (17) 

Then, the decay rate can be obtained from the amplitude (0) as 

rik^l,^) = ^\CLLfml. (18) 

By using the m'j^ ^j given in eq. ([l3|) , we have calculated the branching ratios of the processes 
T ^ fi'j and fi — > 67- We have checked that for mG2 = 45GeV, the branching ratios for 
each processes have the maximums 0(10^^) and 0(10^^^), respectively, at the possible 
minimum value of the slepton mass ~ 45GeV, which are, however, much smaller than the 
present experimental limits (Br(r — > /i7) < 4.2 x 10"^ and Br(/i — 67) < 4.9 x 10^^^ 0]). 

So far, we have neglected the left-right mixings in the slepton mass matrix. However, 
if the SUSY invariant mass fin for Higgs doublets is so big that m^/i// tan/? ~ 0{ml), 
we have non-negligible left-right mixings in the slepton mass matrix. As we will show 
below, these mixing terms give rise to contributions to the Zj Ij'-f process larger than 
the previous estimate in eq.(|18]). The main reason for this is that the chirality flip in the 
fermion line occurs at the internal gaugino mass term (see fig. |^) while in the previous 
case it occurs at the external lepton mass term. Since the gaugino mass is much bigger 
than the lepton mass, these new diagrams may yield dominant contributions. 

Before calculating the diagrams in fig. ^, we first diagonalize the charged slepton mass 
matrix given by the following 6x6 matrix; 

C^.ss = -{L\^)('^k """M f it V (19) 



with 



^Imj = rnl^. + m\cos2^[^m^ew - -jSij, (20) 
^'RR,ij = rn\..-m\cos2[3sm^ew5ij, (21) 
'^Li?,ii = m.^ijiiH taxi (5 + AE,jiV cos (5/ \/2 (22) 



^ We have assumed to; . ^ mi. 
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where mz is the Z-boson mass and sin^ Oy/ the Weinberg angle. Notice that as far as the 
a-parameter in eq.([TT|) is 0(1), the AE^ij is neghgibly small in the case of hh^sji(3 ^ rriQ 
which concerns us. m\ at the electroweak scale are determined by solving the RGEs 



in the same way of the previous calculation for mi^ij, and are found to be m|;_ 
(1.0 - lA)mlSij. 

With the diagonalization matrix U for eq. (|T9D and the eigenvalues (^ = 1^6), 



we can write the amplitude for the process Zj —>■ Ij'-f as 

M{h ^ Ijj) = ClrI^PrU, 4]k, (23) 

where the coefficient Clr is given by 

Clr = T^e^i"^Gi UjAUii+3^-TLR{mlJml^), (24) 

iovr ^ m^j^ 

with 

We find that in the case of small this amplitude ( plj) is well approximated by 

CLRL^,_^^,ertion = ^[^^dl^Gl E ^ LjkUlki{mi^liH tau^ + Ae^uV COS (3 / ^2) 

1 



16vr2 



X 



1 ^ / mL \ 1 ^ / m^i 



— —J-LR — — 2 — ^LR — 2 



em \ em / SLk \ eLk ■ 

2 



(26) 



which is obtained by using the mass insertion. Here, Ul is the diagonalization 

matrix for the mass matrix of the left-handed sleptons (pO]), Tn^^^ its eigenvalues and 
m?^. = m?iif^^^ in eq. (PT|) . From the amplitude (P^D, we get the decay rate as 

T{k^l,^) = ^\CLR\'ml (27) 

which should be compared with the result in eq.(|l8D. 

For a given set of tan/3 and mG2, the decay rates depend basically on the parameters 
and tjiq. However, it is convenient to express them in terms of two physical parameters, 
and {f^if^ < fnf2)j which are mass eigenvalues for the sleptons belonging mainly 

to the r-family. The results are shown in fig. R| for tan j3 = 3, with 17102 = 45GeV and 



7 



rrif-^ = 50GeV being fixed. For a comparison, we also show the results for the case of 
tan P = 30. We easily see that the branching ratios for r — > and /i — 67 are predicted 

as 

Br(r ^ /i7) ~ 3 x 10"^ - 4 x 10^^ (28) 
Br(/i^e7) ~ 5x10^^2-2x10"", (29) 

for 771^2 — (100 ~ 250)GeV.[] It should be stressed that the large branching ratios in 
eqs.([28|) and (^) are never obtained in the case where the left-right mass- insert ion is 
applicable (/i//tan/? ^ mo). To see how the results depend on m^^ and the gaugino mass, 
we show our results in fig. ^ taking rrif^ = lOOGeV and mG2 = 90GeV. From figs. ^ and 
^, we see that the obtained branching ratios are much larger than the previous estimates 
from eq. ([l^) and they lie in the range which will be studied experimentally in no distant 
future. 

We should note that our results shown in figs. ^ and ^ are also larger than the SUSY- 
GUT predictions 0. This is because the authors in ref. have not taken into account 
the left-right mixing effects. However, if these effects are induced, the similar conclusion 
may be obtained even in the SUSY SU(5) GUTs. 

We now briefly discuss other lepton-fiavor violating processes such as /i ^ 3e and 
r ^ 3yU. With the parameter space discussed in the present paper, there also appears 
an enhancement factor mG2/^T in the Penguin diagrams, which will give dominant con- 
tributions to these processes. In this case, the decay rates of /i — > 3e and r ^ 3/i have 
simple relations to those of the fi ^ e-y and r ^ fij processes as 

Br(^-.3e)|pe ^,.„ ^ ^ In ^ ^ 0.8 x 10^, (30) 

Br(r-.3/.)|pe ^ ^ in ^ 0.4 x 10^, (31) 
Br(r — i> /i7) 37r 

where we have taken only the logarithmic contributions. As for the /i-e conversion in 
nuclei, the amplitudes of the box diagrams depend on the squark masses. Thus, in the 
large-squark-mass region, the box diagrams yields negligible contributions. The detailed 



^The minimum value {mf^ ~ lOOGeV) comes from the constraint that the hghtest sneutrino should 
be heavier than 41.7GeV [H. 
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analysis on various lepton-flavor violating processes including the /x-e conversion will be 
given in the future publication |T4 . 



Several comments are in order. 
i) We have considered the parameter space where m^^ ~ (50 — 100)GeV ~ (100 — 
250)GeV and tan/3 = 3 — 30. These parameters correspond to the original parameters 
fiH and mo as 

^ r (l-5)TeV : tan/5 = 3, , . 

\ (100- 500)GeV : tan/5 = 30, ^ ' 

and 

mo = (100 - 250)GeV : tan /3 = 3 - 30. (33) 



With such large ^h-, the radiative electroweak symmetry breaking [1^ hardly occurs as 



long as the universal soft SUSY-breaking mass in eq. (|I0|) is imposed. A solution to this 
difficulty is to abandon the GUT-like relation among the gaugino masses in eq. (|T2|) so 
that the gluino gives rise to much larger masses for squarks than those for sleptons at the 
electroweak scale. 

ii) One may expect that the hierarchy in the stau masses (m^^ <^ ^f2) will give a large 
contribution to p-parameter. We find, however, that their contribution is not substantial 
as far as the mass of the heavier stau is less than 500 GeV for nif^ = 50GeV. 
in) The large fin tan j3 induces also a large left-right mixing in the sbottom sector. We find 
that in some of the parameter space in eqs. ( |32D and (^), the lighter sbottom mass m^^ 
may come in the region excluded experimentally (m^^ ^45GeV). However, this problem 
can be easily solved by giving the gluino a mass larger than expected from the GUT-like 
relation (0), since in this case the lighter sbottom can be lifted above 45GeV through the 
radiative corrections. This is also favorable for suppressing the 6 — > 57 decay substantially 



as will be discussed in ref. 

iv) As pointed out in ref. [Q, the similar lepton-fiavor violation may also occur in the 
framework of SUSY SU(5) GUTs. It should be noted here that in the present model 
the soft SUSY breaking masses m| for left-handed sleptons receive significant radiative 
corrections from the new Yukawa couplings yu,ij, while in the SUSY SU(5) case, those for 
right-handed sleptons, Tn\^p are subject to the large renormalization effects. Thus, in 
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the present model, the i = j = 3 element of m'j^ ^j becomes smaller than other diagonal 
elements, mj^^^ and m|225 by amount of O((10 — 30)%) [1^ as was shown in eq.(0), 
whereas in the SUSY SU(5) GUTs the similar mass shift appears in the right-handed 
slepton sector. Therefore, if the soft SUSY-breaking parameters are precisely determined 
in future experiments, these two scenarios may be distinguished. 
Note added 



After completing this work we became aware that in the very recent paper the 
lepton-fiavor violation due to the Yukawa couplings y,y^ij are also examined in the context 
of SUSY SO(IO) grand unification. The k — > Ij-y process is calculated there by using the 
m|j:j mass-insertion, and our formula in eq. p6| ) is consistent with their result of the y,^^ij 
Yukawa coupling effects. This is, however, only an approximation of our full formula (P3) 
in the small region. 
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Figure 1: Feynman diagrams which give rise to k Ij^y. In each diagram, the blob 
indicates the flavor-violating mass insertion of the left-handed slepton and at the cross 
mark the external lepton flips its chirality. The symbols en, uu, B, W3, and W~ represent 
left-handed charged sleptons, left-handed sneutrinos, bino, neutral wino, and charged 
wino, respectively. 



Figure 2: Feynman diagram which gives rise to Ij'y. The blobs indicate insertions of 
the flavor- violating mass {fn'j^-j) and the left-right mixing mass {fn\iiii)i and at the cross 
mark chirality flip of the bino [B) occurs. The symbols en and em represent left-handed 
charged sleptons and right-handed charged sleptons, respectively. Notice, however, that 
we do not use the mass-insertion method in our calculation as stressed in the text. 



Figure 3: Branching ratios for the processes r — *• /X7 and /i — *■ 67 as functions of mf^. 
The solid lines correspond to i?r(r — > /i7) and the dashed lines to Br{ix ej). Here, 
we have taken mG2 = 45GeV and rrif^ = 50GeV. We also show the present experimental 
upper bounds for each processes by the solid lines with hatches. 



Figure 4: Same as fig. ^ except for 17102 = 90GeV and rrir^ = lOOGeV. 
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